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Abstract 
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1 Introduction 



The atmospheric neutrino (AN) flux is one of the main tools for studying 
neutrino oscillations, neutrino decay and nonstandard neutrino-matter inter- 
actions. Furthermore it represents an unavoidable background for key exper- 
iments in astroparticle physics, as those searching for proton decay, n — n 
transitions in nuclei, etc. Therefore, the accuracy of theoretical predictions of 
the AN flux ultimately affects the interpretation of many experimental data. 

After several studies performed in the frameworks of one-dimensional (ID) [1- 
3] and three-dimensional (3D) [4-6] cascade models, the physics of neutrino 
production in the sub- and multi-GeV energy ranges is fairly well understood. 
However, the validity of the AN flux calculations is still controversial. This is 
mainly related to the uncertainties in the required input data (inclusive cross 
sections for particle production, primary cosmic ray spectrum and composi- 
tion, etc.). Moreover, the problem is too complex to be solved without many 
simplifications (explicit or implicit) and these also are essential sources of un- 
certainties in the predicted AN fluxes [7~10]. A popular opinion is that the 
accuracy of the present-day calculations of the low-energy AN flux is about 
20-30%. Sometimes this belief of obscure origin is treated as the uncertainty 
in the normalization factor. But, as a matter of fact, the uncertainties in the 
input data are so large that the actual accuracy of the AN flux turns out 
to be much worse. It is also essential that the uncertainties are very energy 
dependent and therefore they cannot be reduced to a simple renormalization. 

In this work we present a new ID calculation of low and intermediate energy 
muon and neutrino fluxes, based on up-to-date data on primary cosmic-ray 
flux and hadronic interactions. Our calculation is based on a kinetic approach. 
In order to check the validity of our description of hadronic interactions and 
shower development, we perform a comparison of the predicted atmospheric 
muon fluxes with the data from several recent balloon-borne experiments. 

We believe that an accurate ID AN flux calculation is still helpful, despite the 
results of 3D calculations are now available. Indeed, the main lessons we got 
from the recent 3D analyses [5,6,8,9] may be summarized as follows. 

• 3D effects lead to a strong modification of zenith and azimuth angle dis- 
tributions of neutrinos at energies E^, < 0.5 GeV. On the other hand, the 
corresponding effect for the neutrino-induced events in underground de- 
tectors is expected to be modest and scarcely measurable, because their 
distributions are close to isotropic out of a very big mean scattering angle 
in quasielastic neutrino-nucleus collisions at low energies. 

• Even at very low energies the 3D corrections to the 47r and even 27r aver- 
aged AN fluxes are relatively small. At least, they are comparable with the 
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uncertainties in the input data, including the uncertainties in the cross sec- 
tions for exclusive neutrino-nucleus interactions needed for simulating the 
neutrino-induced events in the underground detectors. 
• Above 1 — 2 GeV the 3D effects become negligible for the neutrino fluxes 
averaged over the azimuth angles at any zenith angle. 

On that ground, we expect the calculation offered for consideration to be valid 
without restrictions at Ei, > 1 GeV and we consider it as a basis for future 
3D generalization of our approach to the AN problem. 

Furthermore, our ID code is useful for testing ID and 3D codes with a method 
completely different from Monte Carlo. It can be used to perform a system- 
atic analysis of uncertainties in the AN flux calculations and to adjust some 
poorly known input parameters, using the data on cosmic-ray secondaries. In 
this work, in order to evaluate quantitatively the uncertainties related to the 
particle interactions, we perform a comparison of the results obtained with 
different models for meson production in nucleon-nucleus and nucleus-nucleus 
collisions. 



2 The CORT code 



The present work is based on an updated code CORT (Cosmic- Origin Radiation 
Transport), previously developed in [11] and used in [1] (see also [12,13]). Like 
the earlier version, the current Fortran 90 code implements a numerical solu- 
tion of a system of ID kinetic equations describing the propagation of nuclei, 
nucleons, light mesons, muons and (anti)neutrinos of low and intermediate 
energies through a spherical, nonisothermal atmosphere. It takes into account 
solar modulation and geomagnetic effects, energy loss of charged particles, 
muon polarization and depolarization effect. The exact kinematics is used in 
description of particle interactions and decays. Let us sketch some most dis- 
tinctive aspects. 

In order to evaluate geomagnetic effects and to take into account the anisotropy 
of the primary cosmic-ray flux in the vicinity of the Earth, we use the method 
of ref. [11] and detailed maps of the effective vertical cutoff rigidities from ref. 
[14] The maps are corrected for the geomagnetic pole drift and compared 
with the later results reviewed by Smart and Shea [15] and with the recent data 
on the proton flux in near earth orbit obtained with the Alpha Magnetic Spec- 
trometer (AMS) [16]. The interpolation between the reference points of the 
maps is performed by means of two-dimensional local B-spline. The Quenby- 
Wenk relation (see e.g. [14]), re-normalized to the vertical cutoffs, is applied 

^ The penumbral effects are included in the definition of the effective cutoff. 
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for evaluating the effective cutoffs for oblique directions. More sophisticated 
effects, like the short-period variations of the geomagnetic field, Forbush de- 
crease, re-entrant cosmic-ray albedo contribution, etc., are neglected. We also 
neglect the geomagnetic bending of the trajectories of charged secondaries 
and multiple scattering effects. Validity of our treatment of secondary nucle- 
ons and nuclei was confirmed using all available data on proton spectra in the 
atmosphere (see [11]). 

The meteorological effects are included using the Dorman model of the atmo- 
sphere [17] which assumes an isothermal stratosphere and constant gradient of 
temperature (as a function of depths) below the tropopause. Ionization, radia- 
tive and photonuclear muon energy losses are treated as continuous processes. 
This approximation is quite tolerable for atmospheric depths h <2x 10^ g/ cm^ 
at all energies of interest [18]. Propagation of /i"*" and fi~ originating from ev- 
ery source (pion or kaon decay) is described by separate kinetic equations for 
muons with definite polarization at production. These equations automatically 
account for muon depolarization through the energy loss (but not through the 
Coulomb scattering). 

In the new version, essentially all approximate formulas used in [1,11-13] are 
replaced (or duplicated for a control) with the more accurate numerical equiv- 
alents. 



3 Primary cosmic ray spectrum and composition 



In the present calculations, the nuclear component of primary cosmic rays is 
broken up into 5 principal groups: H, He, CNO, Ne-S and Fe with average 
atomic masses A of 1, 4, 15, 27 and 56, respectively. We do not take into 
account the isotopic composition of the primary nuclei and assume Z = A/2 
for A > 1, since the expected effect on the secondary lepton fluxes is esti- 
mated to be small with respect to present-day experimental uncertainties in 
the absolute cosmic-ray flux and chemical composition.!^ 

We use the following 9-parameter model for the differential energy spectra of 
every nuclear group at the top of the atmosphere 

dF^iE) ^ I Fief -^'-^ at E < Wa, 



° However, a more rigorous analysis must allow at least for small fractions of pri- 
mary and ^He affecting the /Lt"^//i~ and i^/V ratios. 
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with 

rA(6A) = E it^ ln*= eA, eA - r^E/W^. (2) 

fe=0 

Here £^ is the total energy in GeV/nucleon. The condition for continuity of 
the spectrum and its first derivative in the point E — Wa yields 

In {Fl/Fl) = E ^-7?^ In'^-^^ ta, 7a = E(A; + 1)7?^ In'^ rA- 

k=l k=0 

Thus, we remain with 7 independent parameters for each group which we 
determine by fitting the experimental data. 

For the hydrogen and hehum groups at < 120 GeV/nucleon we chose the 
data of the balloon-borne experiment BESS obtained by a flight in 1998 [19] (a 
period close to a minimum of solar activity). For higher energies (but below the 
knee) we use data by a series of twelve balloon flights of JACEE [20] and the re- 
sult of an analysis by Wiebel-Sooth et al. [21] based upon a very representative 
compilation of world data on primaries and solid theoretical considerations. 
Since the fits obtained in [21] turned out to be very close to the combined 
result of the JACEE experiments, we call our model "BESS-I- JACEE" fit. 
A maximum likelihood analysis with Wi = W4 = 2 TeV/nucleon gives the 
following values for the remaining parameters in eqs. (1) and (2): 
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where are in m~^s^^sr^^(GeV/nucleon)^^. Clearly the many digits in the 
above parameters are by no means indicating the accuracy of the data, but 
they are necessary to avoid loss of accuracy in numerical calculations. 

Figure 1 shows a comparison between the BESS+JACEE fit, the data from 
[19,20], the fit from [21] (shaded areas) and several other experiments [22- 
32], performed in different periods of solar activity. The filled/shaded areas in 
fig. 1 represent the power-type parametrizations derived in[20,22,24,28,31,32] 
from the original data. The legend indicates the publication dates and (when 
relevant) the years of measurements. 

One sees that the BESS-I- JACEE fit is in excellent agreement with the recent 



5 



c 
o 

"2 
s 



10 



10' 



X 



70" - 



10 - 



1 t 



Hydrogen 



Ryan et al., 1972 
Ichimura et al, 1993 
Asakimori et al, 1998 (JACEE) 
Wiebel-Sooth et al, 1998 
Aharonian et al, 1999 (HEGRA) 
Alcaraz et al, 2000 (AMS) 
BESS+JACEE best fit 




ft Seo et al, 1991 (IMPS/LEAP 87) 
A Buckley et al, 1994 (RICH-1 91) 
T Diehl et al, 1997 (RICH-II 96) 

* Menu et al, 1997 (IMAX 92) 
a Bellotti et al, 1999 (MASS 91) 

o Boezio et al, 1999 (CAPRICE 94) 

• Sanuki et al, 2000 (BESS 98) 
a Alcaraz et al, 2000 (AMS 98) 



10 10 



10 

E, ( GeV/nucleon ) 



Fig. 1. Differential proton and helium spectra as a function of kinetic energy per 
nucleon. The data points and filled/shaded areas are from [19-32]. Solid curves 
represent the BESS+JACEE best fit (see text). 

AMS data on the proton flux [31]. On the other hand it shghtly exceeds the 
AMS hehum flux [32] (with average discrepancy of about 15%). Our choice of 
the BESS 98 data was in particular conditioned by the fact that the power-law 
extrapolation of the high-energy tail of the AMS helium rigidity spectrum (cx 
^-(2.74±o.o2)^ [32] distinctly worse joins the world-averaged fit (cx E-(2-64±o.o2)^ 

[21] and even the JACEE 1-12 fit (oc E-i^-^^-o°ot)) [20].g 



We assume that the spectra of the remaining three nuclear groups are similar 
to the hehum spectrum, 

dFA(E) dFJE) 

^ ' CA — (A = 15,27,56). 



dE 



dE 



This is all the more true for the remaining data presented in fig. 1. 
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This assumption does not contradict the world data for the CNO and Ne-S 
nuclear groups but works a bit worse for the iron group. Nevertheless, a more 
sophisticated model would be unpractical since the corresponding correction 
would affect the secondary lepton fluxes by a negligible margin. A normaliza- 
tion using the relevant data yields C15 = 0.068, C27 = 0.026, c^q = 0.0054. 

In this paper we do not consider the effects of solar modulation and use the 
BESS+JACEE spectrum in all calculations. Therefore the predicted muon 
and neutrino fluxes are to some extent the maximum ones possible within our 
ID cascade model. 



4 Nucleon-nucleus and nucleus-nucleus interactions 

Calculations with the earlier version of CORT [1,11-13] were based on semiem- 
pirical models for inclusive nucleon and light meson production in collisions of 
nucleons with nuclei, proposed by Kimel' and Mokhov (KM) [33] and by Serov 
and Sychev (SS) [34] (see also [35,36] for the most recent versions). The KM 
model is valid for projectile nucleon momenta above ~ 4 GeV/c and for the 
secondary nucleon, pion and kaon momenta above 450, 150 and 300 MeV/c, 
respectively. Outside these ranges (that is mainly within the region of reso- 
nance production of pions) the SS model was used. Comparison of the KM 
model with the data and with some other models is discussed in [7,37]. 

Both models are in essence comprehensive parametrizations of the relevant 
accelerator data. In our opinion, the combined "KM+SS" model provides a 
rather safe and model-independent basis to the low-energy atmospheric muon 
and neutrino calculations. However it is not free of uncertainties. For the 
present calculation, the fitting parameters of the KM model for meson and 
nucleon production off different nuclear targets were updated using acceler- 
ator data not available for the original analysis [33,35]. The values of the 
parameters were extrapolated to the air nuclei (N, O, Ar, C).\[\ The overall 
correction is less than 10-15% within the kinematic regions significant to at- 
mospheric cascade calculations. Besides, energy-dependent correction factors 
were introduced into the model to tune up the output /n^ ratio taking into 
account the relevant new data. This correction only slightly affects the total 
meson multiplicities and modifies the pion charge ratio within 5-9%. 

The processes of meson regeneration and charge exchange (vr'^-l-Air — n^^^^ + 
X etc.) are not of critical importance for production of leptons with energies 
of our interest and can be considered in a simplified way. Here we use a proper 
renormalization of the meson interaction lengths, which was deduced from the 

^ The old version of CORT used the renormalized A^Be inclusive cross sections. 
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results of ref. [38] obtained for high-energy cascades. 



The next important ingredient of any cascade calculations is a model for 
nucleus-nucleus collisions. The overwhelming majority of muon and neutrino 
flux predictions is based on the simplest superposition model (SM) in which 
the collision of a nucleus A, with total energy and atomic weight A, against 
a target nucleus B, is treated as the superposition of A independent collisions of 
nucleons with the target nucleus, each nucleon having an energy Eq = E^/A. 
This approximation is based on the hypothesis that, when the energy per 
nucleon of the projectile is much larger than the single nucleon binding en- 
ergy, the A nucleons interact incoherently. Generally speaking, this approach 
is rather far from reality and its applicability to the lepton flux calculations is 
not evident. A more accurate treatment of nucleus- nucleus collisions becomes 
especially important at low energies and for the regions or directions with high 
geomagnetic cutoffs. Indeed, the magnetic rigidity of a proton bounded in a 
nucleus is a factor A/Z larger than that of a free proton of the same energy and 
therefore (considering the steepness of the primary spectrum) mainly nuclei 
are responsible for production of the low-energy secondaries. 

Here we consider a modest generalization of a simple "Glauber-like" model 
used in [1,11,12]. Namely, we write the inclusive spectrum of secondary parti- 
cles c (c = p, n, n"^, , K^, . . . ) produced in AB collisions as 



dN. 



—SAB 



dx I dx dx 

+ (1 - Cab) [Z5cp + iA- Z)5,n] 5{l - x). (3) 

Here dNj^-B^cx/ dx is the spectrum of particles c produced in A^B collisions 
(A^ = p, n) and Cab is the average fraction of inelastically interacting nucleons 
of the projectile nucleus A. The term proportional to delta function describes 
the contribution of "spectator" nucleons from the projectile nucleus. 

In the standard Glauber-Gribov multiple scattering theory the quantity 
is certainly independent of the type of inclusive particle c. On the other hand, 
it depends of the type of nucleus collision. Indeed, essentially all nucleons 
participate in the central AB collisions {Cab — 1)Q while, according to the 
well-known Bialas-Bleszyhski-Czyz (BBC) relation [39], 

/•c „inel / „incl f a\ 

^AB = O^AfB/^AB (4) 

for the minimum bias collisions. 



To use (3) in a cascade calculation one should take into account that nu- 
cleons and mesons are effectively produced in nuclear collisions of different 

^ Here we suppose for simplicity that the atomic weight of the projectile nucleus is 
not much larger than that of the target nucleus. 
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kind. Namely, the contribution from central collisions is almost inessential for 
the nucleon component of the cascade but quite important for light meson 
production. Thus one can expect that effectively ^X'b^ > ^^b- 

From now on we use relation (4) for nucleon production by any nucleus while 
for meson production we put ^He-Air = ^ ^^^Y ^^e parameter ^ from 0.517 
(the BBC value) to 0.710 (an experimental upper limit derived from the data 
on interactions of a particles with light nuclei [40]). Similar (within ±15%) 
variations of the parameters ^X'^ir for other nuclear groups lead to very small 
effects on the muon and neutrino fluxes; we do not discuss these effects below. 

It is evident that the SM corresponds to the formal limit ^^b = 1 in eq. (3), if 
one simultaneously puts cr^g = cr]^3 in the transport equations for the nuclei 
and mesons. Below, we use this approximation making efforts to represent the 
results of the Bartol code [2] which is based on the TARGET model for meson 
production and the SM for collisions of nuclei. 



5 Numerical results: muon fluxes 

In order to test the validity of our description of hadronic interactions and 
shower development, we compare in this section some numerical results for 
muon fluxes with data from several recent balloon-borne experiments. More- 
over, in order to evaluate quantitatively theoretical uncertainties, we compare 
the calculations obtained by using different codes and/or interaction models. 
Namely, we consider numerical results obtained with: 
i) the CORT code, using our "standard" (KM+SS) interaction model; 
a) the CORT code, using the TARGET model for tt/K meson production^ 
(below, we refer to this calculation as "CORT+TARGET"). 
Hi) the original Bartol Monte Carlo code [2] with the same (BESS+JACEE) 
fit for the primary spectra ("TARGET"). 

In fig. 2 we compare the fi~ fluxes calculated with the three mentioned models 
for twelve atmospheric ranges with the data of the balloon-borne experiment 
CAPRICE 94 [41]. The measurements were performed during ground runs in 
Lynn Lake {h = 10^ g/cm^) and during the continuous ascension of the balloon 
to its float altitude {h ^ 4 g/cm^). The nominal geomagnetic cutoff rigidity 
Rc was about 0.5 GV and the detection cone opening was of about 20° around 
the vertical direction with the typical zenith angle of 9°. The values of depths 
indicated in fig. 2 are the Flux-weighted Average Depths (FAD) [41] for eleven 

^ Due to technical reasons it is difficult to adopt also the TARGET model for 
nucleon production; we thus describe the nucleon production by using the KM 
model everywhere. 
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atmospheric ranges A/ij {hi < 10^ g/cm^).P^ 




Muon Momentum ( GeV/c ) Muon Momentum ( GeV/c ) 

Fig. 2. Differential momentum spectra of negative muons for different atmospheric 
depths. The data points are from the balloon-borne experiment CAPRICE 94 [41]. 
The curves and filled areas are the results of calculations with CORT and TAR- 
GET (see text). The numbers over the curves indicate the flux- weighted average 
atmospheric depths (in g/cm^) and scale factors. 

Figure 2 contains many data. In order to make it easily readable, let us explain 
the notation in detail. The solid curves represent the results of our calculations 
for each FAD. They are obtained with the standard CORT and with ^ varying 
between the limits mentioned in section 4, that is within ±15%. The thickness 
of the curves reflects the indetermination in the parameter ^. For obvious 
reasons, the muon flux uncertainty due to this indetermination is maximal 

There are also small (less than 3%) spreads of the FADs within some of the ranges 
Ahi, but these are completely negligible. 
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at the top of the atmosphere (it is about 15% at h = 3.9 g/cm^ and p = 
100 MeV/c) and becomes almost neghgible for h > 100 g/cm^. The wide 
filled areas display the variations of the muon fluxes inside the ranges A/ij. 
It is important to note that the thickness of the bands is relatively small 
just for the region of effective muon and neutrino production that is in the 
neighborhood of the broad maxima of the muon flux (100 — 300 g/cm^). This 
means that, in this region, the evaluation of the FAD cannot introduce any 
relevant uncertainty. Outside the region of effective production of leptons, 
the amplitude of the muon flux variations increases with decreasing muon 
momenta on account for the strong dependence of the meson production rate 
upon the depth aX h < 100 g/cm^ and the growing role of the muon energy loss 
and decay at /i > 300 g/cm^. The dashed curves show the results obtained 
using the original Bartol code and the dotted curves are the results of the 
CORT+TARGET model (we shall comment them later). 

Figure 3 shows a comparison of the calculated differential momentum spec- 
tra of yU+ and n~ with the ground level data obtained in the experiments 
CAPRICE 94 [41,42] {h = 1000 g/cm^, = 0.5 GV) and CAPRICE 97 [42] 
{h — 886 g/cm^, Rc — 4.2 GV). The detection cone in both experiments was 
the same as described above. The calculation are done using the KM+SS in- 
teraction model. Variation of C, leads to a small (< 3%) effect. One sees quite 
a good agreement between the theory and data. 

For comparison, the result obtained using the SM is also shown. The data be- 
low ~ 10 GeV/c are precise enough in order to conclude that they disfavor the 
superposition model. Let us note that the calculated sea-level {h = 10^ g/cm^) 
/x"*" -|- /x~ flux is systematically lower than that obtained using the old ver- 
sion of CORT (see ref. [13]): the difference is about 18% for p = 10 GeV/c 
and vanishes only at p > 80 GeV/c. The main reason of this difference is 
that the BESS-I-JACEE primary spectrum is essentially lower than the spec- 
trum used in ref. [13] and in the earlier calculations [1,11] for energies below 
~ 200 GeV/nucleon. Another reasons pertain to the changes in the inclusive 
cross sections and improved description of muon propagation. 

In figs. 4 and 5 we present the muon charge ratio fi~^ / fi~ as a function of muon 
momentum for different atmospheric depths and geomagnetic cutoffs. The 
data are from several balloon experiments: MASS 89 [43], MASS 91 [48,47], 
IMAX92 [44], HEAT 94 [49], HEAT 95 (two data handlings of the same data 
sample given in [45] and [46]), CAPRICE 94 [41] and CAPRICE 97 [42]. The 
world average best fit, recently obtained in ref. [50] by an analysis of the data 
from many ground level experiments at p > 10 GeV/c, is also shown. The 
muon charge ratio, being sensitive to the primary chemical composition and 
to the tt^/tt" ratio at production, is a useful tool for testing the predicted v/u 
ratio. 
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10 1 10 10 

p (GeV/c) 

Fig. 3. Differential momentum spectra of positive and negative muons at the 
ground level. The data points are from the experiments CAPRICE 94 [41,42] and 
CAPRICE 97 [42]. The curves are calculated with CORT. 
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Fig. 4. Muon charge ratio near the top of the atmosphere. The points and shaded 
rectangles represent the data from the experiments MASS 91 [47], 1MAX92 [44], 
HEAT 95 [45,46] and CAPRICE 94 [41]. The curves and filled areas are the results 
of calculations with CORT and TARGET. 



Several important conclusions can be deduced by comparing the various data 
presented in figs. 2, 3, 4 and 5. 

• There is a substantial agreement between the "standard" CORT predictions 
and the current muon data within wide ranges of muon momenta and atmo- 
spheric depths. In particular, the agreement is good for the region of effective 
production of leptons, in which the spread of the data (partially related to 
the procedure of determination of the FAD and mean muon momenta) is 
minimal. This provides an evidence for the validity of our model. 

• The comparison between CORT and CORT+TARGET allows to elucidate 
the "hadronic" uncertainties in muon fluxes. The only difference between 
the two calculations is indeed the assumed meson production model and 
the model for nucleus-nucleus interactions. The ii~ fluxes obtained using 
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Fig. 5. Muon charge ratio at the ground level. The points and shaded rectangles 
represent the data from the experiments MASS 89 [43], MASS 91 [48], IMAX92 [44], 
HEAT 94 [49], HEAT 95 [46], CAPRICE 94 [41] and CAPRICE 97 [42]. The curves 
and filled areas represent the results of calculations with CORT and TARGET and 
the best fit obtained in ref. [50] from many experiments at p > 10 GeV/c. 



the TARGET model systematically exceed the results obtained with our 
standard model. In order to have a quantitative estimate of the hadronic 
uncertainties, we note that the muon flux at p = 1 GeV/c and h = (100 — 
300) g/cm^ (i.e. the range where muons and neutrinos are effectively pro- 
duced) increases by about 50% when the TARGET model is appUed. For 
comparison, we remind that the accuracy of the CAPRICE 94 experiment 
in this region is about 15%. 
• We can see that, except for very small depths, there is a serious disagree- 
ment between the CORT+TARGET and the original Bartol results. The 
disagreement grows fast with decreasing muon momentum and with in- 
creasing depth, and it is comparable (or even larger) to the hadronic uncer- 
tainties. Since the two calculations are obtained by using the same primary 
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spectrum and the same meson production model, the observed difference 
can only be due to different features in the nuclear cascade development and 
muon propagation. On the other hand, the observed agreement at floating 
altitude, where such features are not important, assures the correctness of 
our use of the TARGET meson production model. Possible sources for the 
observed difference could be nucleon elasticity distributions, meson regen- 
eration processes, muon energy losses and some more tiny features of the 
cascade models under consideration. 



6 Numerical results: neutrino fluxes 

In this section we present our numerical results for atmospheric neutrino (AN) 
flux. Due to geomagnetic effects, the AN spectra and angular distributions are 
quite different for different underground neutrino experiments. However, for 
the aims of present study it is enough to consider only one representative case. 
Let us confine ourselves to results pertinent to the Kamioka site. 

Figure 6 shows the z/g, l/g, z/^ and 17^ energy spectra averaged over all zenith 
and azimuth angles. The shaded areas are the results obtained with CORT 
using our standard interaction model. The widths of the areas indicate the un- 
certainty due to variations of the ^ parameter. One sees that this uncertainty 
is at most 6% and thus it is negligible. The dashed curves correspond to the 
CORT+TARGET model, while the circles show the results by a new (prelim- 
inary as yet) 3D calculation based on the FLUKA Monte Carlo simulation 
package [5] (see also [9]). Note that the primary spectrum used in the latter 
calculation is also a parametrization of the BESS 98 data, which is however 
not identical to our BESS-fJACEE fit. In fig. 7 we compare the zenith-angle 
distributions of z/g, ^e, and V^, calculated with CORT (adopting ^ = 0.685), 
CORT+TARGET and FLUKA; the distributions are averaged over azimuth 
angles and over seven energy ranges. The comparison allows to "highlight" 
the 3D effects which are highly dependent on neutrino energy and direction 
of arrival. 

As a first point, we note that our standard model leads to the lowest AN 
flux at low energies. Our flux is also essentially lower than that predicted by 
Bartol group [2] and by Honda et al. [3]p^ and that are used in many analy- 
ses of the data from 1MB, Frejus, NUSEX, Kamiokande, Super-Kamiokande, 
SOUDAN 2 and MACRO. 

We do not enter here in the complex problem of the interpretation of the AN 

Let us remind that the results of refs. [2,3] were obtained with different (higher) 
fluxes of the primary cosmic rays. 
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Fig. 6. Scaled in averaged fluxes of I'ejT^e, z^/x and 17^ for Kamioka site. Shaded areas 
represent the result of CORT obtained with ^ varying between 0.517 and 0.710. The 
dashed curves are the result of CORT obtained with the TARGET model for meson 
production and superposition model for collisions of nuclei. The circles are for the 
result of a 3D calculation based on the FLUKA code [5] . 



anomaly in the light of our result. We want to remark, however, that, in the 
context of our model, it is difficult to increase the AN flux without spoiling the 
agreement with the current data on hadronic interactions, primary spectrum 
and muon fluxes. Let us emphasize that our low AN flux, when it is applied to 
the analysis of the underground neutrino experiments, results in some electron 
excess together with (or rather than) the muon deflcit in the neutrino induced 
events. 

Let us now discuss in detail the differences between our results and those 
obtained within CORT+TARGET and FLUKA. The quantitative comparison 
is given in tables 1 and 2. Table 1 shows the ratios of in averaged AN fluxes. 
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Fig. 7. Zenith angle distributions of i^e, i^ei t'/i and 17^ for several energy ranges. 
Solid curves represent the result of CORT obtained with = 0.685 and the rest 
notation is the same as in fig. 6. 

obtained with CORT+TARGET and FLUKA to those with CORT for several 
energy ranges. In table 2 we tabulate the "flavor ratio" , defined by 
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Table 1 

47r averaged AN fluxes for Kamioka site calculated with CORT+TARGET and 
FLUKA and normalized to CORT. 





CORT+TARGET 




FLUKA 




(GeV) 


















0.1-0.2 


1.64 


1.78 


1.72 


1.73 


1.25 


1.45 


1.41 


1.40 


0.2-0.3 


1.51 


1.68 


1.60 


1.61 


1.23 


1.38 


1.34 


1.33 


0.3-0.5 


1.42 


1.61 


1.48 


1.50 


1.17 


1.28 


1.22 


1.22 


0.5-0.7 


1.34 


1.50 


1.36 


1.38 


1.10 


1.19 


1.11 


1.12 


0.7-1.0 


1.28 


1.41 


1.27 


1.30 


1.05 


1.12 


1.04 


1.05 


1.0-2.0 


1.20 


1.30 


1.17 


1.20 


1.02 


1.05 


0.99 


1.01 


2.0-3.0 


1.10 


1.17 


1.03 


1.08 


0.98 


1.03 


0.96 


0.97 


3.0-5.0 


1.02 


1.07 


0.95 


1.01 


0.99 


1.02 


0.95 


0.97 


5.0-7.0 


0.95 


0.99 


0.91 


0.96 


1.01 


1.05 


0.94 


0.99 


7.0-10. 


0.91 


0.95 


0.89 


0.93 


1.03 


1.09 


0.95 


1.04 


10.-20. 


0.87 


0.91 


0.88 


0.91 


1.07 


1.12 


0.96 


1.02 



where i/g, Fe, etc. stand for the At: averaged fluxes. This quantity is represen- 
tative for the ratio of e like to like single-ring contained events measured in 
water Cherenkov detectors and to "showers-to-tracks" ratio measured in iron 
detectors. 

The comparison between CORT and CORT -h TARGET predictions provides 
a direct determination of the "hadronic" uncertainties in the AN fluxes. We 
remind indeed that these two results are obtained by using the same primary 
fluxes and the same code; thus the observed differences is only due to the 
assumed meson production model. At low energies, the CORT+TARGET 
predictions drastically exceed those from CORT and even at energies about 
(2-2.5) GeV one can observe a sizeable difference between the two results. At 
~ 0.5 GeV (tab. 1), the discrepancy is of the order of 40-50%. This shows 
that the hadronic uncertainties are quite relevant in the energy range which 
is responsible for the sub-GeV event rate in the Super-Kamiokande detector. 

Figs. 6 and 7 show that CORT and FLUKA neutrino fluxes are in good agree- 
ment for any zenith angle at energies above ~ 1 GeV. On the other hand, 
CORT is systematically lower than FLUKA for < 1 GeV. At the energy 
Ej^ = 0.5 GeV the FLUKA An averaged z/g, Vf., i^fj, and 17^ fluxes exceed the 
corresponding CORT results by about 15%, 27%, 18% and 18%, respectively. 
This discrepancy is only partially due to 3D effects (fig. 7), which can account 
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Table 2 

Neutrino flavor ratios R^, calculated with CORT, CORT+TARGET and FLUKA 
for Kamioka site. 



A^;^ (GeV) 


CORT 


CORT+TARGET 


FLUKA 


0.1-0.2 


0.52 


0.51 


0.48 


0.2-0.3 


0.52 


0.50 


0.49 


0.3-0.5 


0.51 


0.50 


0.50 


0.5-0.7 


0.50 


0.50 


0.50 


0.7-1.0 


0.49 


0.50 


0.50 


1.0-2.0 


0.47 


0.49 


0.49 


2.0-3.0 


0.44 


0.47 


0.45 


3.0-5.0 


0.40 


0.43 


0.42 


5.0-7.0 


0.36 


0.37 


0.38 


7.0-10. 


0.32 


0.32 


0.34 


10.-20. 


0.26 


0.26 


0.29 



for an increase < 10%, and it is most probably related to differences between 
the FLUKA and CORT hadronic interaction models. 

Finally, as it is seen from table 2, the flavor ratio is essentially the same in all 
three models, except for very low energies. As one can conclude from fig. 7, 
the same is also true for the "up-to-down" asymmetry, that is for the ratio of 
upward-going (cos ^ < 0) and downward-going (cos ^ > 0) neutrino fluxes. 



7 Conclusions 

Let us summarize here the main points of our discussion. 

• Calculations with a new version of the ID code CORT, using the updated 
KM+SS interaction model and the BESS+JACEE primary spectrum are 
in agreement with the recent data on muon momentum spectra and charge 
ratios measured at different atmospheric depths and geomagnetic locations. 

• The uncertainty due to indetermination of the parameters of our model for 
nucleus-nucleus interactions is comparatively small everywhere, except for 
high altitudes. Some experimental uncertainties (like those due to determi- 
nation of the flux-weighted average depths) are negligible for the regions 
of effective generation of leptons and at ground level. Therefore the agree- 
ment between the predictions and the muon data obtained at these depths 
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provides a conclusive evidence for validity of our approach. 

• Our atmospheric neutrino flux is systematically lower than those used in 
the current analyses of the data from underground experiments, while the 
neutrino flavor ratio and "up-to-down" asymmetry are essentially the same 
as in all recent calculations. We stress that, in the context of our model, it 
is difficult to increase the AN flux without spoiling the agreement with the 
current data on hadronic interactions, primary spectrum and muon fluxes. 

In order to estimate theoretical uncertainties in the muon and neutrino fluxes 
we compared our results with those obtained by using different codes and/or 
interaction models, with the following conclusions. 

• By implementing two different interaction models within the same code, we 
have shown that the "hadronic" uncertainties are quite large at low energies 
for both muon and neutrino fluxes. Specifically, by using the TARGET 
meson production model instead of the KM+SS, we obtained a sizeable 
increase of neutrino fluxes at energies up to ~ 3 GeV. At the energy ~ 
0.5 GeV representative for the sub-GeV events in the Super-Kamiokande 
detector, one observes a 40-50% increase of the AN fluxes. 

• By comparing two different codes implementing the same interaction model 
(TARGET), we have shown that other sources of theoretical uncertainties 
are also relevant. Namely, a sizeable difference exists between the muon 
fluxes predicted by the original Bartol code and those obtained by imple- 
menting the TARGET interaction model within the CORT code. We believe 
that the observed disagreement is due to the differences in the calculation 
of nuclear cascade, muon energy loss and decay. 

• A comparison of the 47r averaged AN fluxes calculated with CORT and 
FLUKA shows rather good agreement at Ei, > 0.7 — 0.8 GeV (the difference 
is typically less than 15%) and a sizable disagreement at lower energies. A 
comparison of neutrino zenith-angle distributions, predicted by CORT and 
FLUKA, suggests that the main sources of the disagreement at low energies 
are the differences in the adopted interaction models and 3D effects. 

Finally we conclude that the AN spectra and angular distributions predicted 
with CORT using the BESS+JACEE primary spectrum and KM-SS interac- 
tion model are quite safe at energies above ~ 1 GeV while the low-energy 
range (which is of special importance, as a source of background to the proton 
decay experiments) is still uncertain. 
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